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ABSTRACT 

The Vulpecula OB association, VulOBl, is a region of active star formation located in the 
Galactic plane at 2.3 kpc from the Sun. Previous studies suggest that sequential star forma- 
tion is propagating along this 100 pc long molecular complex. In this paper, we use Spitzer 
MIPSGAL and GLIMPSE data to reconstruct the star formation history of VulOBl, and search 
for signatures of past triggering events. We make a census of Young Stellar Objects (YSO) in 
VulOBl based on IR color and magnitude criteria, and we rely on the properties and nature of 
these YSOs to trace recent episodes of massive star formation. We find 856 YSO candidates, 
and show that the evolutionary stage of the YSO population in VulOBl is rather homogeneous 
- ruling out the scenario of propagating star formation. We estimate the current star formation 
efficiency to be ^ 8%. We also report the discovery of a dozen pillar-like structures, which are 
confirmed to be sites of small scale triggered star formation. 

Subject headings: Infrared: ISM, Stars — ISM: HII regions — Stars: Formation, Pre-main sequence 



1. Introduction 

The physical mechanisms describing stellar 
birth are fairly well understood for low- and 
intermediate-mass stars but still under debate 



for th eir high-mass analogues (|McKee &: Ostrikei 
20071 ). On larger scales, observations show that 
star forming mechanisms are of relatively poor ef- 
ficiency, as only a fraction of the gas reservoir in 
the Universe is turning into stars. Typical Star 
Formation Efficienci es (SEE) are of th e order 3- 
6% in the Galaxy (lEvans et al.l l2009l). a nd 5% 
or less in other galaxies ( Rownd fc Yound 19991) . 
Still, in extreme environments like in the star- 
burst galaxy Arp 220, the S EE can reach 50% 
( Anantharamaiah et al. I I2OOOI) suggesting that the 
star forming mechanism at work in such objects 
could be of a different nature. Eor instance the 
feedback into the interstellar medium (ISM) from 
short-lived massive stars seems to influence the 



yield of star formation in their local environment. 
As they evolve off the main sequence, high-mass 
stars produce a copious amount of energy while 
still embedded in their native cocoon; such dis- 
ruption of a molecular cloud leads to gravitational 
instabilities and possibly t o the onset of a new 
episod e of star formation. iHosokawa &: Inutsuka 
(j2006l ) showed that under certain conditions run- 
away triggering can take place around massive 
OB stars. In other cases however, turbulence and 
magnetic fields can have a negative feedback on 
the local ISM and lea d to the suppression of the 
star fo rming activ i tv (iPrice fc Batd [200! IStond 
19701) . lElmegreenI (Il998l) provides an exhaustive 
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review of the theoretical framework to study the 
dynamical triggers of star formation. 

In this context, we use th e data from the Sp itzer 
Legacy surveys MIPSGAL ('Carcv et al1l2009[ ) and 
GLIMPSE (,Beniamin ct al. 200a) to study the 
star forming activity currently taking place in 
the Vulpecula OB association (hereafter Vul OBI). 
VulOBl hosts nearly one hundred OB stars and 
three bright Hll regi ons known as S harpless ob- 
jects 86, 87 and 88 (jSharplessI [l959l) . According 
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to lEhlerova et alj pOOll ) andlTurner| (|l986l ). the 
star forming activity occurring in VulOBl might 
have been triggered by a common external source, 
and star formation might be propagating from one 
Hll region to another through the expansion of a 
supernova shock front. The aim of this paper is 
to study the triggered star formation on scales as 
large as ~100 pc in VulOBl, which could help 
us understand star formation mechanisms within 
other Giant Molecular Clouds in the Galaxy and 
beyond. 



We use the method developed bv lGutermuth et al 
io obtain a census of Young Stellar Ob- 
jects (hereafter YSOs) in VulOBl based on in- 
frared color and magnitude criteria, and we rely 
on the properties and nature of these YSOs to 
trace episodes of star formation. The MIPSGAL 
and GLIMPSE sensitivity makes our search for 
YSOs in VulOBl biased towards massive objects 
that reveal the most recent star forming activity. 

In Section [2] we present the dataset used in the 
paper. In Section |3] we give a comprehensive de- 
scription of VulOBl encompassing the three Hll 
regions as well as the dozen pillar-like structures 
we discovered in this region. Section 0] describes 
the identification process of Young Stellar Objects 
based on their IR-excess emission. We present and 
discuss our results in Section [SJ and we give our 
conclusions in Section |6l 



2.1. 



The dataset 

Spitzer Observations and Point Source 
Catalogs 



The Spitzer Space Telescope ([Werner et al 



20041 ) observed over 270 square degrees of the inner 
Galactic plane in six wavelength bands as part of 
two legacy programs: the Galactic Legacy Infrared 
Mid-Plane Survey Extraordinaire (GLIMPSE; 



E. Churchwell PI; iBeniamin erahl [2003) and 



the MIPS G ALactic plane su rvey (MIPSGAL; 
S. Carey PI; ICarev et alll2009l) . The Vulpecula 
region was covered by GLIMPSE I (PID 188) 
and MIPSGAL I (PID 20597). It was imaged 
with the IRAC camera (jPazio et all 12004 ) at 3.6, 
4.5, 5.8 and 8.0 M m, and with the MIPS camera 
(jRieke et al.ll2004l) at 24 and 70 fim. The angular 
resolution of Spitzer is 2", 6" and 18" at 8, 24 and 
70 /xm, respectively. 



was used to build mosaics of the sky about 
1 square degree wide from individual Basic Cali- 
brated Data (BCD) frames. Details of the post- 
processing a nd BCD pipe l ine n iodifications are 
describ ed in iMeade et al.l (l2007l) for GLIMPSE 



and in iMizuno et al.l ()2008f ) for MIPSGAL. The 



large Vulpecula mosaics presented in this paper 
are also built with MOPEX using the 1 square de- 
gree plates available through the Infrared Science 
ArchiveQ (IRS A). 

The search for YSO candidates in Vulpecula 
relies on the point source catalogs (PSC) gener- 
ated from the GLIMPSE and MIPSGAL surveys. 
For GLIMPSE, point sources are extracted using 
a modified version of the Point- S pread Functio n 
(PSF) fitting program DAOPHOT (|Stetsonlll987l) . 
The GLIMPSE PSC we used is the enhanced 
data product v2.0 available on the IRSA web- 
site. For MIPSGAL, sources are extracted from 
the 24 /im maps using a Pixel- Response Function 
fitting method (Shenoy et al. in preparation) and 
the APEX software developed at the Spitzer Sci- 
ence Cente'^. This catalog is then merged with 
the GLIMPSE PSC to include IRAC and 2MASdl 
fluxes if they are available. 70 fim sources are ex- 
tracted from reprocessed 70 fim maps (Paladini et 
al. in preparation) u sing the StarFinder software 
(|Diolaiti et al.ll2000l ). and matched to the MIPS- 
GAL PSC with a search radius of 9". 

2.2. Ancillary data 

We complement our dataset with visible, sub- 
millimeter, millimeter and radio data to obtain a 
broader view of the Vulpecula complex, and to in- 
vestigate the morphology, star formation history 
and possible relationship between the Hll regions. 
We make use of the VLA Galactic Plane Survey 
(VGPS) data, which consists of Hi line and 21- 
cm continuum emission data with a resolution of 



1^x1 ^x1.56 km s ^ and 1' respectively (iStil et al 



20061 ). These data are analysed in details in Sec- 



The MOPEX package (jMakovoz et al.l I2006D 



lhttp://irsa.ipac.caltech.edu/data/SPITZER/GLIMPSE/ 
and 

htt p://irsa.ipac.caltech.edu/data/SPITZER/MIPSGAL7] 
^http://ssc.spitzer. caltech.edu/postbcd/apex.html 
^The Two Micron All Sky Survey (2MASS) is a joint project 
of the University of Massachusetts and the Infrared Pro- 
cessing and Analysis Center/California Institute of Tech- 
nology, funded by the National Aeronautics and Space Ad- 
ministration and the National Science Foundation. 
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tion 13.11 We also use the Virginia Tech Spectral- 
line Survey (VTSS) which provides arcminute- 
resolution images of the 6563 A Ha recombina- 
tion li ne of atomic Hydroge n in the Vulpecula re- 
gion (jPennison et al.l Il998 ). S. Bontemps (pri- 



Schneider et al.l 



r s (pri- 
(l2006h 



vate communication, see 
for details) also provided our team with a 1.2'- 
resolution extinction map of VulOBl generated 
from stellar counts and colors derived from the 
2MASS point source catalog (see Figure [J). The 
visual extinction in Vulpecula ranges from 3 to 
20 with a distribution peaking at Av~5. We also 
use the CO surve y of the Galaxy published in 
(|20nih . 



Dame et al 



However the relatively low 
spatial resolution of ~7.5' over Vul OBI prevented 
a relevant detailed study of the CO emission. We 
find nonetheless a very good correlation between 
the CO emission and the visual extinction that 
both trace the cold and dusty molecular mate- 
rial (cf Figure [1]) . In addition, we used the sub- 
millimeter data obtained during the second flight 
of the Balloon-borne L arge Aperture Subm illime- 
ter Telescope (BLAST. IPascale et al.lliooih . where 
a 4 square degree region was mapped around Sh2 - 



86 at 250, 350 and 500 /xm (Cha pin et al]l2008[ ) 



The angular resolution of the BLAST maps are 
40", 50" and 60" at 250, 350 and 500 ^im, re- 
spectively. And finally, we exploit the recently re- 
leased 1.1mm data of the Galactic Plane Survey 
carried out with Caltech Submillimeter Observa- 
tory/Bolocam (Aguirre et al. in preparation) and 
providing an angular resolution of 31". 

3. The Vulpecula OB association 

VulOBl is an active star forming region lo- 
cated in the Galactic plane at a longitude l'^60°. 
In this article we focus on the region centered at 
(Z, b) — (60.2, —0.2) and cover ing about 6.6 square 
degree s (see Figures [Hand [2]). ICarmanv fc Stencel 
(|1992[ ) describe VulOBl as an oval region about 
3.5° by 1.5° hosting the star cluster NGC6823. 
Nearly one hundred hot massive stars are found 
in the direction of Vulp ecula f rom t he catalog 
of OB stars compiled by iReedl (|2003h . Some of 
these stars have already been associated with and 
are responsible for the ionization of three bright 
Hll regions in Vul OBl, namely Sh2-86, 87 and 
88 (|Sharplessiri959t ). Figure [2] presents the mid- 
infrared morphology of VulOBl and the location 
of the OB stars. 



The VGPS 21-cm continuum emission cov- 
ering VulOBl reveals two bright circular com- 
pact sources (^^1.3' in diameter) coincident with 
Sh2-87 and Sh2-88, and an extended dimmer re- 
gion (27'xl2') reminiscent of the infrared mor- 
phology of Sh2 -86. The super nova remnant 



SNR G59.5-h0.1 (|Tavlor et al.lll992l ) is also clearly 
visible in the radio continuum data at [l, b) = 
(59.59,0.1), its measured diameter is 15' (cf Sec- 
tion 15.4.11) . Our analysis of the VGPS Hi line 
data indicates that the three Hll regions are 
neighbors (cf Se c tion 13.11), as was pre v iously 

(|200lh 



not ed by Turner ( 198(1) . Ehlerova et al 



and lCappa et al.l (|2002n . Distance determinations 



ing stars range from 1.5 to 3.2 kpc (Fich & Blitz 




1984 Barsonvl 1989 


; Guetteij 1992; Brand & Blitz 




I993I iMassev et al. 


1995: Deharveng et al. 2000l 




Hovle et al. 2003: Kharchenko et al.ll2005:lBica et al. 


2008h. Following: the arguments fromlChaoin et al. 



(20081), we adopt a common distance of 2.3 kpc 
for the three Hll regions. 

In the catalog of iReedl (120031) we find about 



30 OB stars in the direction of Sh2-86 and Sh2- 
87, in what we identify as the pillars region (these 
stars might be responsible for the existence of 
pillar-structures, see Section 1X2)) . We check their 
possible association with VulOBl by computing 
their distance moduli based on the My-Spectral 
type relation and the i ntrinsic colors given by 



Landolt-BornsteinI ( 1982f) . We assume a ratio of 



total-to-selective absorption i?v — 3.1 typical of 
the diffuse ISM, and we associate a visual extinc- 
tion to each star based on our extinction map 
of VulOBl. Out of the 16 stars for which spec- 
tral type and m agnitude information are avail- 



able in the IReedl catalog, 12 fall between 1.3 and 
2.3 kpc with a mean distance from the Sun of 
1.91 ± 0.4 kpc. We find similar distances using 
intrinsic colors and a bsolute niagnit udes given by 
Johnson! (|l958l ) and iWegneil (|2006l ) respectively. 



Although the uncertainties are significant, which 
is mostly due to the large uncertainty on the values 
of the absolute magnitudes, our distance estimate 
for these stars is consistent with an association 
with the Vulpecula complex. 

3.1. Hll regions 

Sh2-86 is an ex tended Hll region about 40' wide 
(ISharplesd Il959l) located in the southern part 
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Fig. 1. — CO contours derived from the velocity- integrated CO map ofE ame et al.l (j200l[ ). integrated from 
20 to 40 km s~^, overlaid on a visual extinction map of VulOBl (Av ranges from 3 to 20 mag on a linear 
scale). The Sharpless Objects appear as higher extinction regions. 



of Vulpecula (Figures [U [2] and [TT] in the ap- 
pendices). It is excited by the open cluster 
NGC 682 3, which is part of t he OB associatio n 



VulOBl dMassev et al.l Il995l: iBica et al.l |2008[) 



Pigulski et all (|2000[ ) estimate the age of the star 



cluster at 3 ± 1 Myr old. In the atomic gas dis- 
tribution derived from the VGPS data, Sh2-86 
appears as an oblong hole at velocities between 
26 and 31 km s~^. We find that the dimension 
of the shell is approximately 60' along the Galac- 
tic plane and 25' across it. In the optical, the 
south-eastern part of Sh2-86 is very contrasted. It 
exhibits silhouetted pillar-like structures pointing 
towards NGC 6823, emission and reflection nebu- 
lae as well as filamentary structures. The infrared 
emission is also contrasted in this region tracing 
the complex interplay between the UV radiation 
and the surrounding dusty interstellar medium. 
Chapin et al.l ( 20081 ) also report the detection of 



49 compact sub- millimeter sources associated with 
Sh2-86 using the Balloon-borne Large Aperture 
Submillimeter Telescope (BLAST). The presence 
of these clumps, with masses ranging from 14 to 
700 Mq , indicates the capability of Sh2-86 to form 



massive stars. 

Sh2-87 an d Sh2-88 are also Hll regions initially 
discovered bv lSharpless ( 19591 ). These regions are 
active sites of star formation as indicated by the 
presence of H2 O maser hne emissions and bipola r 
outflows ( Barsonv 19891 : Deharveng et al.l l2000f l . 
Bolocam data show bright sources at the loca- 
tion of Sh2-87 and Sh2-88 indicating the pres- 
ence of cold core candidates. Our analysis of the 
VGPS data reveals that the morphology of the Hi 
shell around Sh2-87 at velocities between 22 and 
26 km is very similar to that of the 24 ^m emis- 
sion. The Hll region appears as a hole in the Hi 
distribution as t he hydrogen is ionize d by an em- 
bedded BO star (|Felh fc Hartenlll98ll ). A compact 
source seen in absorption sits at the center of Sh2- 
87, its position is coincident with the location of 
the peak emissi on at 24 urn (I, b) = (60 .88, —0.13). 
Barsonvl (|l989f l and IXue fc Wul (l2008h give a de- 



tailed description of this Hll region. 

Sh2-88 is a diffuse nebula of diameter ~20' lo- 
cated at {l,b) = (61.45, 0.34), which is exc i ted b y 



the 08 star BD -f25°3952 (|Cappa et al.ll2002h . 



The star formation activity taking place in Sh2- 
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88 occurs in a couple o f neb ular knots identified 
bv iLortet-Zuckermannl (|l974l ). Sh2-88A and Sh2- 
88B, located 15' southeast of Sli2-88. These knots 
are responsible for the bright 24 /j,m emission vis- 
ible in Figure Hat {l,b) = (61.47,0.1). Sh2-88B 
is actually the brightest of the three Sharpless ob- 
jects at mid-IR wavelengths. Figure [10] in the ap- 
pendices shows a composite image of Sh-88B in 
IRAC bands. It consists of a compact cometary 
Hll region and a n ultra compact (UC) Hll region. 
Deharveng et al. I (I2OOOI) present a detailed analy- 
sis of Sh2-88B and its stellar content. At radio 
wavelength, the association of Sh2-88 with a Hi 
shell is somewhat more difficult than in the other 
two cases as the edges of the shell are not as well 
defined. We detect however a faint hole at veloci- 
ties between 22 and 26 km as well as a compact 
source seen in absorption in the center of Sh2-88B. 

3.2. Pillar structures 

We have discovered several pillar-like struc- 
tures on either side of the Galactic equator be- 
tween Sh2-86 and Sh2-87 (59°.5 < / < 61°). 
These objects are similar to the archetypal pil- 
lars of creation found in M16 (iHester et al.lll996l: 
Urguhart et al. I l2003l) . Such pillars are usually as- 
sociated with recent episodes of star formation as 
the winds and radiation emanating from young 
massive stars are responsible for sculpting these 
elongated elephant trunks out of the surrounding 
molecular material. An accepted mechanism for 
the formation of the pillars is the slow photoevap- 
oration of a pre-existing molecular clump shad- 
owing a tail of more diffuse gas, but other mech- 
anisms have been proposed base d on hydrody - 
namical i nstabi lit ies for instance. Spitzen (|l954 ) , 
Bertoldil (Il989[) . iLefloch fc Lazarefj (|l994[) and 
Carlqvist et al. ( 20031 ) present analvtical mod- 
els for the for r nation of pillar-like struct ures, 
and iMiao et al.l (l2006l). iMizuta et al.l (|2006l ) and 



Gritschneder et al. 



( 2009t ) carry out numerical 
simulations of the formation and evolution of such 
objects. According to these studies, the pressure 
at the surface of the pillars due to the strong ex- 
ternal radiation appears to trigger the formation 
of new stars, which is confirmed by recent observa 

ri u^^j ^ ■■ ' 



tions fe.g. ISug itani et al. 2002; 



Re ach et al, 
.2001) 



2004 : 



ture to identiiy individual objects. The three pil- 
lars VulP12-13-14 were already know n prior to 
these observations ( Chapin et al.l I2OO8I ) : they are 
located northeast of NCG 6823 and are obviously 
associated with NGC 6823 as they point towards 
the star cluster. However the 11 pillars VulPl 
to VulPll have never been reported in the lit- 
erature before this study. Their association with 
VulOBl is not as straightforward as in the case 
of NGC 6823 since no obvious source could be 
identified as their sculptor. The OB stars men- 
tioned previously are certainly good candidates, 
for instance the star HD 186746 (spectral type 
B8Ia) is located right above the tip of the pillar 
VulPl in projection (Figures [5] and [3]) . Neverthe- 
less, 10 of the discovered pillars, out of 11, seem 
to point towards the same faint diffuse nebulos- 
ity located at {I, h) — (60.37, —0.04), which is not 
coincident with any known OB stars. We argue 
that the chance for ten randomly distributed pil- 
lars to point towards the same object is negligi- 
blfl thus the discovered pillars VulPl - VulPl are 
most certainly associated with each other. The 
central nebulosity might host the object respon- 
sible for the formation of the pillars. In addi- 
tion, a diffuse faint circularly symmetric structure 
centered on the same nebulosity is discernible in 
the background diffuse emission. Bright objects 
at 24 /im also seem to delineate this faint struc- 
ture at a radius of 15' centered at (^, 6) = 
(60.37, —0.04). These morphological features sug- 
gest that a single event in the past may have 
molded the pillars region. 

We analyze the radio data in the pillars region 
to look for a possible association between these 
structures and the OB association. Out of the 
15 pillars identified in Figure[3l only 3 are detected 
in the VGPS Hi hue data. The pillar VulPlO 
presents a deep absorption feature at 27-33 km s~^ 
as well as smaller features centered at 8, 13 and 
20 km s~^. These Hi features are marginally re- 
solved into two close compact sources located at 
the base of the pillar. Since these features are 
quite different from the IR morphology of VulPlO, 
we cannot exclude a fortuitous association with a 
coincident radio source seen in projection. 

VulP4 and VulP5 appear to have similar char- 



Bowler et al.ll2009; .Reach et al.. 

Figure |3] presents the 24 fj,m image of the pil- 
lars found in Vulpecula, and gives the nomencla- 



^Assuming we measure the direction of the pillars with an 
accuracy of ~ 2°, the probability is lower than 10^^^. 
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Fig. 2. — The top panel shows a composite image of VulOBl (3.6, 8.0 and 24 fj,m are color-coded as blue, 
green and red respectively). This image was assembled by R. Hurt and T. Robitaille based on GLIMPSE and 
MIPSGAL data. The bottom panel represents Ha contours from the VTSS overlaid on the MIPSGAL 24 ^m 
image. The Ha emission associated with Sh2-88 peaks at {l,b) — (61.45,0.30) while the bright mid-IR 
emission located 15' south-east of the diffuse nebula is the knot identified in the text as Sh2-88B. Green 
crosses indicate the location of the 88 OB stars present in this 3.7° x 1.8° region. 
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acteristics, they present comparable morphologies 
in the Hi line data and in the Spitzer bands. They 
both emerge from the molecular cloud to be ex- 
posed to the ionizing radiation. These pillars are 
located at the edge of Hi shells, and their silhou- 
ette is clearly identified at velocities 21-23 km s~^, 
which places them at roughly the same distance 
as the three Sharpless objects for which Vlsr = 
22 -31 km s-\ 

VGPS continuum emission shows a bright two- 
lobed compact source located at the base of the 
pillar VulPlO, consistent with the two blobs ob- 
served in the Hi line data. No other pillars are 
detected in the 21-cm continuum. 

Tables [3] and H] in the appendices summarize 
the geometrical and morphological information ex- 
tracted from the Spitzer observations for individ- 
ual pillars. Of particular interest are the bright 
(~ 2 Jy) compact sources found at 70 /im in the 
core of the pillars VulPl and VulPS (see Fig- 
ure |4]). These embedded sources are faint at 
24 /xm 1 mjy) and have no counterparts in 
IRAC and 2MASS bands. Even though star for- 
mation is expected to occur preferentially at the 
tip of such structures, these red sources are likely 
sites of massive star formation. Longer wave- 
lengtlil observations would be necessary to iden- 
tify the nature of these sources. Figure|4]shows the 
pillars VulPS and VulP6 and their associated red 
sources. YSO candidates are located at the tip 
of each pillar, and also in the pedestal of VulPS 
along some Infrared Dark Clouds (IRDC). Sub- 
millimeter BLAST sources are also located at the 
tip of the pillars and along the same IRDCs. Two 
Bolocam sources are also located along one of the 
IRDCs. 

4. The IR-excess emission of Young Stel- 
lar Objects 

Stars form from the gravitational collapse and 
fragmentation of giant molecular clouds in the 
ISM. The contraction of the cold gas leads to 
the formation of a dense rotating core radiating 
its thermal energy in the millimeter and far-IR 
regime. As the bulk of the initial cloud mass 
is falling towards the center of the core, a pro- 



^The BLAST data do not cover these pillars, and the 1.1 mm 
Bolocam data, which is available for these pillars, is not as 
good a diagnostic as the sub-millimeter data. 



tostar emerges with a gaseous and dusty accre- 
tion disk rotating around it. The peak emis- 
sion of such a young stellar object shifts towards 
shorter wavelengths revealing a bright source in 
the mid-IR. The circumstellar disk is then dissi- 
pated via accretion, planet formation or evapora- 
tion. Further reading on star formation mecha- 



nisms c an be found i n e.g. [ Terebev et al 
Adams et all (fl987l): lAndre et all (|l993l . 



1984) 



20001) : 



McKee fc Ostrikeil (|2007l) 

Recent works have shown that the wavelength 
coverage of Spitzer instruments is well suited 
for observing the mid-IR excess emission em- 
anating from the warm circumstell ar material 
around Young Stell a r Obj ects 



Guieu et all l2009l) 



20041: iHaryev et all l2007t iKoenig et al 



Allen et al. 
2008t 



B a sed on the mid-IR ex- 
tinction law (iLutzi 11999 : Indebetouw et al. 1 120051: 



Flaherty! l2007l ). ICutermuth et al.l (|2008l ) further- 
more argue that the most reliable IRAC-based 
criterion for identifying YSO candidates is the 
[4.5]-[5.8] color. Indeed the fiattening of the ex- 
tinction curve observed between the IRAC bands 
2 and 4 (4.5 to 8 /im) reduces the degeneracy be- 
tween selective interstellar extinction and intrinsic 
IR excess. 

Our approach to identify YS Os in Vul OBI is to 
exploi t the method developed bv lGutermuth et aT 



( 20081 ) based on Spitzer colors and magnitude cuts 
to identify sources with IR-excess emission. We 
apply this method to the Point Source Catalogs 
(PSC) compiled by the GLIMPSE and MIPS- 
GAL teams (Section 12. 1|) to take full advantage 
of the good photometric quality and the high re- 
liability of these data. The direct consequence 
of this approach is that no extended sources will 
be considered in our study even though they ex- 
hibit significant IR-excess emission. To mitigate 
this caveat, we compare the angular resolution of 
Spitzer with the typical angular size of circum- 
stellar envelops and disks seen at 2.3 kpc from the 
Sun. Assuming the telescope optic is diffraction 
limitecjfl down to 5.5 /im, the resolution element at 
3.6 and 24 /im corresponds to 3500 and 16000 AU 
at 2.3 kpc respectively. The spatial extend of the 
Young Stellar Objects we are loo king for does not 
exceed these limits. For instance, Vicente fc Alved 



^From the Spitzer Observ er's Manual available at 
|http: / /ssc.spitzer.caltech.edu| 
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Fig. 3. — Identification of pillar structures in VulOBl on a MIPSGAL 24 /im image. The grey scale is 
stretched to increase the contrast of the pillars. VulPl to VulPlO seem to point toward the same direction. 
Extended concentric structures are noticeable around the central position up to a radius of 10 pc. VulP12 
to VulP14 are associated with the young star cluster NGC 6823. 



(|2005l ) measured the size distribution of circum- 
stellar disks in the Trapezium cluster and found 
disk sizes ranging from 100 to 400 AU, an order 
of magnitude smaller than the resolution element 
at 3.6 ^m. Although protostellar envelops have 
a large r extend, i.e. a few thous and AU accord- 
i ng to Bonnell et al.l ( 19961) and J0rgensen et al 



(|2002n . and assuming they can only be detected 
at long wavelengths where the resolution element 
is larger, we expect the number of resolved YSOs 
in VulOBl to be relatively small. However a frac- 
tion of the YSOs in their early phases of evolution 
might be resolved by Spitzer, even at 24 /im, and 
those large protostars would be excluded from our 
analysis. 

4.1. Initial catalog 

We first compile a catalog of IRAC point 
sources based on the highly reliable Point Source 
Catalog available from the IRSA website as a 
GLIMPSE I v2.0 enhanced data produclQ. All 
sources found in this catalog with galactic coordi- 
nates in the range 58.6 ° <l < 62.2 ° and -1.0 ° < 



& < 0.8° are considered 3.28 x lO'^ sources). 
Fluxes are converted to magnitudes using the ze- 
ropoints given in Table [TJ To ensure good pho- 
tometry and avoid flux stealin g betwe e n adj acent 



sources as noted by Robitaille et al.l ()2008D . we 



require that valid point sources have magnitude 
errors lower than 0.2 mag, and that their close 
source flail (csf) is zero. We also apply specific 
cuts on IRAC magnitudes to achieve high detec- 
tion reliability, the magnitude limits we adopt are 
14.2, 14.1, 11.9 and 10.8 mag at 3.6, 4.5, 5.8 and 
8.0 /im respectively. According to the GLIMPSE I 
Assurance Quality Documeni|£|, these limits ensure 
a detection reliability above 98% in all four IRAC 
bands. The reliability of the faintest sources how- 
ever might vary across the VulOBl region as it 
depends on the background brightness structure 
and the local source density. 

This intermediate catalog is then merged 
with the MIPSGAL 24 /im Point Source Cata- 



""http: / /www. astro, wise. cdu/glimpse/glimpsel_dataprod_v2.0.pdf 



Following the eonvention of the GLIMPSE I v2.0 enhanced 
data products, csf=0 implies that no sources in the Archive 
Catalog are within 3' of the source, 
^http: / /www. astro. wisc.edu/glimpse/GQA-master.pdf 
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Fig. 4. — Spitzer 3-color images of some pillars discovered in VulOBl. YSO candidates are marked as white 
small squares and circles, large blue boxes show BLAST sources and green ellipses indicate Bolocam sources. 
Top panel: Five pillars identified as VulPll and VulPl to VulP4 in Figure [3] (from left to right). Blue, 
green and red are 8, 24 and 70 /im images, respectively. VulPl and VulP2 exhibit bright 70 /im sources in 
their core. Bottom panel: VulP5 and VulP6. Blue, green and red are respectively 4.5, 8.0 and 24 /im images. 
InfraRed Dark Clouds (IRDCs) are visible in the pedestal of the pillar VulPS. YSO candidates are aligned 
along IRDCs. 
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log (Shenoy et al. in preparation) with the con- 
straint that the maximum separation between two 
matching sources is 3". Once again to ensure a 
good photometry, we require that the magnitude 
of 24 /im sources and their associated errors are 
constrained. These constraints however are ap- 
phed at a later stage of the selection procedure 
(see Section 112]) • 

Out of nearly 3.3 x 10^ sources present in the 
GLIMPSE catalog within our search area, only ^ 
1.47 X 10^ bypassed the initial magnitude cuts. 
Most of the exclusions were due to noisy detections 
in bands 3 and 4 of IRAC. 

In the end, the initial catalog from which we 
search for YSO candidates contains point sources 
with well measured photometry and with at least 
one detection in IRAC bands, plus the associated 
2MASS and MIPSGAL 24 /im fluxes if they exist. 

4.2. Mid-IR selection of YSO candidates 

We start our search for YSO candidates in the 
initial catalog by removing extragalactic contam- 
inants that might be misidentified as YSOs. Us- 
ing the selection criteria presen ted in the appen- 
dices of [Cutermuth et al. ( 2008 ). which are based 
on a compl ementary analysis of the Bootes field 
IRAC data (|Stern et al.ll2005[ ). we identify and re- 
ject 6 sources dominated by Polycyclic Aromatic 
Hydrocarbons (PAH) feature emission, likely star- 
forming galaxies or weak-line active galactic nuclei 
(AGN). No sources matching the colors and mag- 
nitudes criteria for broad-line AGNs were found. 
Considering a search area of 6.5 square degrees, we 
find a density of extragalactic sources in Vulpecula 
slightly lower than one per square degree. This 
is comparable to the density of ~0.5 galaxy per 
square degrees found in the zone of avoidance by 
Marleau et aP (l2008l) using GLIMPSE and MIPS- 
G AL data. We furth er filter the catalog according 
to iGutermuth et al.l criteria and reject 8 sources 
that have a large 4.5 /im excess emission consistent 
with molecular hydrogen line emission found in 
regions where high- velocity outflows interact with 
the cold molecular cloud. The remaining sources 
are presumably of stellar origin. 

All sources with the following color constraints 
are considered likely YSOs: 

[4.5] - [8.0] > 0.5, 

[3.6] - [5.8] > 0.35, 



[3.6] - [5.8] < 3.5 X ([4.5] - [8.0]) - 1.25, 

Out of the ~ 2.38 x 10"* sources possessing all four 
IRAC magnitudes, we identified 820 YSO candi- 
dates based on their IR-excess emission. We fur- 
ther classify the selected objects according to their 
infrared spect ral index arp = d log{XF\) / d log{X) 
as defined bv iLada ( 1987 ). We specifically com- 
pute the spectral index amAC as the slope of 
the Spectral Energy Distribution (SED) measured 
from 3.6 to 8.0 /xm. Objects with aiuAC > —0.3 
have a flat-ish or rising SED indicating the pres- 
ence of a cold dusty envelop infalling onto a central 
protostar, these are designated class O/I YSOs. 
Objects with —0.3 > ajRAC > —1-6 are classi- 
fied as class II YSOs, these are pre-main-sequence 
stars with warm optically thick dusty disks or- 
biting around them. We define class HI objects 
as having —1.6 > ajj^AC > —2.56: these stars 
have cleared most of their circumstellar environ- 
ment, their near-IR emission is mostly photo- 
spheric but they may present some excess emis- 
sion above 2 /xm. Such o bject s are dubbed ane- 
mic disks bv iLada et al. ( 2006 ). Finally, objects 
with aiRAc < —2.56 are stars with photospheric 
emission only, this spectral index corresponds to 
the slope of the stellar photosphere SED in the 
Rayleigh- Jeans domain. 

Note that this YSO classification i s based on the 
classif ication scheme proposed by ICreene et al 
(|l994l ). except that the Flat Class (0.3 > a > 



—0.3) is included into the class O/I which repre- 
sents the population of proto stars with infalling 
envelopes. I Gal vet et al.l (|l994l ) indeed showed that 
Flat Class YSOs could be interpreted as infalling 
envelopes. Figure [5] presents a color-color diagram 
based on IRAC data only, the different symbols 
indicate the classes associated with each YSO can- 
didate identified in our analysis. 

The next step is to exploit the extra informa- 
tion contained at longer wavelengths to confirm, 
or reclassify, the YSO candidates identified solely 
from their IRAC colors. We first check that the 
SEDs of all YSO candidates continue to rise from 8 
to 24 /im. We also look for objects previously clas- 
sified as photospheric sources that exhibit bright 
24 /im fluxes ([5.8] - [24] > 1.5). These sources are 
thought to be transition disks, i.e. class II YSOs 
with s ignificant dust clearing from their inne r 
disks (|Calvet et al.l 120021 : iD'Alessio et~al1 12005[ ). 
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Table 1 

Zeropoints used for Flux-to-Magnitude Conversions. 





2MASS 




IRAC 




MIPS 


Band 


J H K 


[3.6] 


[4.5] [5.8] 


[8.0] 


[24] [70] 


Zeropoint [Jy] 


1594 1024 666.7 


280.9 


179.7 115.0 


64.13 


7.17 0.778 



Note.— 2MASS zeropoints are from lCohen et all (j2003l ). IRAC from lReach et all (|2005l ) and 
MIPS from lRieke et al.l (|2008l ). 



these objects are of particular importance for un- 
derstanding the evolution of circumstellar disks 
around young stars. We also check for protostar 
miscla ssifications due to e xtrem e visual extinction 
levels. Gutermuth et al. ( 20081 ) suggests that if a 
protostar that has MIPS detection does not meet 
the criterion [5.8] - [24] > 4 (if they possess [5.8] 
photometry) or [4.5] — [24] > 4, then it is likely a 
highly reddened class II YSO. Finally, we require 
that any sources lacking detections in some IRAC 
bands yet having bright 24 /im fluxes ([24] < 7 and 
[X] — [24] > 4.5 mag, where [X] is the photometry 
for any IRAC detection available in our catalog) 
have to be added to the list of likely YSOs, and to 
be classified as highly embedded protostars. We 
apply the 24^m-based color constraints on YSO 
candidates only if their MIPS detection reliability 
is > 95 % (S. Shenoy, private communication); 
thus we require that, for the MIPS photometry to 
be relevant for the YSO selection, the 24 /im mag- 
nitude is [24] < 8.6 mag and that its associated 
error is (T[24] < 0.2 mag. 

We tested the ability of the method at find- 
ing genuine YSOs on a well studied star forni- 



20051: iHarvev et al.l 12007 
Chavarn'a et al.ll2008l )" 



Robitaille et"al1 l2008t 



ing region, RCW79 (jZavagno et al.l 120061 ). and 



it proved to be adequately reliable and efficient 
(Zavagno, private communication). Furthermore, 
this method has been successfully applied to 
study embedded stella r clusters in NGC 1333 
(jGutermuth et al.ll2008l ) , the star formation activ- 
ity in the Hll region W5 ( Koenig et al.''2008') or in 
the g iant molecular cloud G216-2.5 ([Me geath et aTl 
20091 ). Note that other methods, also based on 
Spitzer colors and magnitudes, have been devel- 
oped to find YSO candidates (e.g. lHartmann et al 



4.3. Final census of YSO candidates 

We finally find 856 YSO candidates in Vul OBI: 
239 are likely protostars with infalling envelops 
(class O/I), 464 are disk-bearing stars (class II), 
and 153 are class III objects with very little cir- 
cumstellar material. Among the class I objects, 
15 YSO candidates are likely deeply embedded 
protostars, and 89 have a spectral index consistent 
with the Flat Class. Among the class II objects, 
85 are highly reddened class II YSOs that were 
misclassified as protostars based on IRAC criteria 
only. A further 21 sources are classified as likely 
transition disk class II objects. 

Except for a few cases, all YSO candidates have 
detections in all four IRAC bands, about half of 
them have 2MASS detections in all three JHK 
bands, 75% have 24 /im counterparts (it reaches 
100% for class O/I objects) and 8% have 70 /im de- 
tections. Table [5] gives the coordinates and fiuxes 
of all YSO candidates. Table [5] in the appendices 
gives the coordinates, fiuxes and type of the con- 
taminants identified in Vul OBI. 

To test the robustness of our analysis relative to 
the interstellar extinction, we have used the visual 
extincti on map of Vul O BI and the formulae pro- 
vided bv lFlahertvl (|2007t ) to deredden IRAC fluxes. 
We run the exact same procedure as described 
above, and we find that the YSO census is only 
marginally different from the uncorrected case (ex- 
cept for the class III population that is reduced by 
35%). This justifies the choice of [5.8]- [8.0] as an 
extinction-free indicator for finding YSOs. 
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Fig. 5. — Left panel: IRAC color-color diagram of point sources found in the Vulpecula region. Stars with 
photospheric emission are marked as black dots in the figure, they have null colors in IRAC bands so they 
aggregate around point (0,0). Color- and symbol-coded YSO candidates protrude from the bulk of the stellar 
sources towards the positiv e color quadran t due to their IR-excess emission. The reddening vector is derived 
from the extinction law of iFlahertv ( 2007 ). Right panel: Color-Magnitude diagram of point sources based 
on MIPS and IRAC bands. (A color version of this figure is available in the online journal.) 



Table 2 

List of YSO candidates and their 2MASS, IRAC and MIPS photometry. 



Galactic 2MASS IRAC MIPS 

GLIMPSE source name Glon Glat J H Ks [3.6] [4.5] [5.8] [8.0] [24] [70] Class 



SSTGLMC G060.3261-00.6407 60.3261 -0.6407 8.68 8.19 7.74 7.07 4.55 ... II 

SSTGLMC G060.1815-00.5632 60.1816 -0.5632 12.71 10.61 9.45 8.31 8.22 7.89 7.42 7.16 ... Ill 

SSTGLMC G060.3304-00.6042 60.3304 -0.6042 12.99 11.40 10.71 10.14 4.78 -5.65 O/I 

SSTGLMC G060.3331-00.6065 60.3331 -0.6065 16.04 14.38 13.41 11.79 10.83 10.12 9.09 5.89 ... O/I 

SSTGLMC G060.2117-00.6700 60.2118 -0.6700 14.69 13.61 13.05 12.30 11.62 11.24 9.95 II 

SSTGLMC G060.2898-00.6427 60.2898 -0.6427 13.14 12.17 11.44 10.74 6.98 ... O/I 

SSTGLMC G060.2868-00.6351 60.2868 -0.6351 12.90 11.92 11.17 10.32 7.65 ... O/I 

SSTGLMC G060.1319-00.5355 60.1320 -0.5355 10.84 8.99 8.01 7.30 7.24 6.94 6.70 5.44 ... Ill 

SSTGLMC G060.2180-00.6362 60.2180 -0.6362 10.68 8.79 7.86 7.24 7.28 6.85 6.55 5.27 ... Ill 

SSTGLMC G060.1975-00.6052 60.1975 -0.6052 13.63 12.32 11.28 10.38 7.24 ... O/I 



Note. — Galactic coordinates are in units of degrees (°). Table [2] is published in its entirety in the electronics edition of the 
Astrophysical Journal. A portion is shown here for guidance regarding its form and content. 
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4.4. Reliability and completeness 

In the present study, we make the deliberate 
choice to favor the highest possible detection re- 
liability, at the expense of a modest completeness 
figure. This choice was largely motivated by the 
need to automate the search for YSOs over such a 
large area of the sky. 

As mentioned in Section 14. 1[ we achieve a de- 
tection reliability better than 98% for the set of 
IRAC magnitude cuts chosen for populating our 
initial catalog. However the actual detection re- 
liability depends to a certain extent on the level 
of the background diffuse emission. An adequate 
assessment of the completeness of our YSOs cat- 
alog would require to carry out a large spectro- 
scopic survey of the red objects in VulOBl in 
order to identify genuine YSOs, and then com- 
pare the results to our catalog. This is impractical 
though, given the large number of (faint) objects 
in VulOBl. 

Considering the stringent criteria used to reach 
a reliability over 98%, we expect a noticeable 
number of YSOs to be excluded from our initial 
point source catalog. For instance, any saturated 
sources, that are either bright sources on faint 
background or faint sources on bright background, 
are not considered in our analysis. There are very 
few saturated data in IRAC images so that the 
YSOs selection is basically unaffected by this ef- 
fect. Still, we expect that a fraction of YSO candi- 
dates will not have 24 /xm counterparts because of 
saturation issues (especially around Sh2-88B and 
Sh2-87 where the background is high). Extended 
sources are also excluded based on the selection 
procedure used in this work (see Section S]), but 
their actual number is expected to be fairly small 
at 2.3 kpc. Point sources with too large a magni- 
tude error are not considered either. We set the 
limit for bad photometry to be Umag > 0.2 mag. 
IRAC bands 3 and 4 being the noisiest bands, the 
selection criterion on CTmag mostly impacts [5.6] 
and [8.0] detections. 

The completeness of our YSO catalog is also 
dependent on the evolutionary stage of the ob- 
ject. For instance, the wavelength band 2-24 /im is 
well suited for detecting and identifying disks 
around young stars (class II), but longer wave- 
lengths are more appropriate for observing early- 
stage objects (class O/I). We believe that our sam- 



ple of IRAC-selected YSO candidates may lack a 
significant fraction of these protostars. For in- 
stance, we looked for a possible correlation be- 
tween the spectral index aiRAC and the visual 
extinction of the YSO candidates expecting to 
find a population of protostars with high values 
of aiRAC: the youngest objects, to be the 

most extinct/embedded objects. However we did 
not find such correlation in our sample which sug- 
gests that we might be missing the most embed- 
ded phases of star formation based on IRAC and 
MIPSGAL 24 fiiR data only. Similarly, the popu- 
lation of class HI YSOs which lacks near-IR excess 
emission is strongly underestimated in our study, 
even if they are genuine YSOs possessing Ha or X- 
ray emission. For all these reasons, it is very diffi- 
cult to assess the completeness of our YSO catalog 
wi th confidence. We u se the evolutionary models 
of iBaraffe et al. ( 19981 ) to compute the mass of a 
2 Myr old star at 2.3 kpc with no IR-excess emis- 
sion having a 3.6 /zm magnitude of 14 mag; and 
we find a mass of ~1 Mq. Considering that the 
8 /um channel of IRAC is the less sensitive chan- 
nel (limit magnitude of 10.8 mag), and that most 
YSOs have significant IR-excess emission, we esti- 
mate our YSO catalog to be complete down to a 
few solar masses. 

Lastly, we have to account for the possible 
misidentification of YSO candidates with evolved 
stars. The circumstellar environment of Plane- 
tary Nebulae (PNe) or Asymptotic Giant Branch 
(AGB) stars are rich in warm dust grains emitting 
in the near- mid-IR such that these objects occupy 
the same locus in color-color diagrams as YSOs 
(|Hora et al.l l2008l: 'Srinivasan ct al."2009l). Follow- 
ing the analysis of Robitaillc et al. (2008!), we es- 
timate the fraction of evolved stars in our cata- 
log of YSO candidates to be about 25%, of which 
most are AGB stars. We will argue in Section [521 
that AGB stars could be further excluded from 
the sample of YSO candidates based on clustering 
criteria. 

5. Results and discussion 
5.1. Star Formation Efficiency 

We estimate the current star formation effi- 
ciency in Vulpecula by comparing the mass of 
the gas reservoir, Mdoud, with the mass that has 
turned into stars during the last few million years. 
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Myso- The Star Formation Efficiency (SFE) is 
then derived as follows: 



SFE: 



Myso 



Myso + Mdoud 



(1) 



We compute the mass of the cloud using our ex- 
tinction map of Vulpecula and the formula relat- 
ing the column density and the visual extinction 
Nh/Av = 1-37 X 10^^ cm~^.mag~^ assuming a 
value of Ry = 5.5 typica l for m olecular clouds as 
suggested in lEvans et al. We use a contour 

of Ay = 7 mag to delineate the active star forming 
region on the Ay map of Vulpecula. Wc find an 
average mass surface density of 16.2 Mq.pc"^. We 
integrate the surface density over the 2.7 x 10'^ pc^ 
enclosing the Ay > 7 region, and we find a cloud 
mass Mcioud = 4.5 x 10^ Mg. 

Since our YSO sample represents only the high- 
end of the Initial Mass Function (IMF), as we es- 
timated in Section we need to account for the 
missing mass of the total YSO population in or- 
der to derive a relevant Myso- We use the IMF 



of Kroupal (j2001l ) to compute the number of stars 
with masses ranging from 0.01 to 50 Mq. We nor- 
malize the IMF assuming 510 stars have masses 
>1 Mq within the Ay > 7 region, and we find 
that VulOBl should contain 2 x 10"* YSOs with a 
total mass of Myso — 4200 Mq. We thus obtain a 
SFE of ~8%, similar to othe r star formin g regions 
in nearby molecular clouds ([Evans et al.t,2009. ). 



5.2. 



Spatial distribution of YSO candi- 
dates 



The spatial distribution of YSO candidates in 
VulOBl is presented in Figure [HI YSO candi- 
dates are not distributed randomly in the field. We 
rather find several coherent structures, or groups 
of YSO candidates, which in most cases are rem- 
iniscent of the mid-infrared morphology of the 
complex. For instance the highest densities of 
YSO candidates are located close to, or within 
the three Hll regions. Other YSOs appear to line 
up along InfraRed Dark Clouds (IRDC), e.g. at 
{l,b) = (59.98,0.06) at the pedestal of the pillar 
VulPS (see Figure S]) ; or along bright contrasted 
structures of the extended 24 fun emission, e.g. at 
{l,b) = (61.82,0.33) or {l,b) = (59.80,0.64). YSO 
candidates have also been identified at the tip of 
most newly discovered pillar structures (see Ta- 
ble H and Sections Eil] and [O:^ . 



If we assume the IMF parameters presented 
in Section 15.11 we estimate that the average sur- 
face density of YSO candidates is 7.4 YSO.pc"^. 
This value is consistent with the typical surface 
density found in other sta r forming regions, e.g. 
13.0 YS0.pc"2 Serpen s (iHarvev et alj|2007l) or 
3.3 YSO.pc-2 in Lupus (|Merin et al.ll2008h . We 
build a surface density map of YSO candidate s 
following the definition of IChavarria et aL (2008). 
At each point of a 5" grid, we compute the surface 
density of YSOs as 

TV 



where is the distance to the = 5 nearest 
neighbor. This image is then convolved with a 
gaussian of FWHM = 1' to obtain a smooth con- 
tour map of the surface density. Figure [5] shows 
the surface density contours over a 70 fiia im- 
age of VulOBl. It clearly delineates the mid- 
IR-bright HlI regions, plus other groups of YSO 
candidates associated with fainter compact mid- 
IR sources, e.g. at (Z,6) = (60.60,-0.70). The 
surface density peaks at ^50 YSO.pc"^ at the 
position {I, b) — (59.80, 0.064) which is coinci- 
dent with the BLAST source V30, also known as 
IRAS 19410-1-2336. This source is a candidate Hy 



per Co mpact HlI region according to lChapin et al 



(j2008[ ). and about 15 of our YSO candidates ap- 
pear to be grouped around it (see Section r5. 4. 2p . 

In addition to this population of clustered YSO 
candidates, the re is a population of dis tributed IR- 



excess sources. 



Robitaille et"aLl (|2008l ) argue that 



a fraction of these isolated red objects could be 
AGB stars mistaken for YSO candidates in our se- 
lection process (they occupy the sa me locus in the 



C-C d iagram, cf Section[4]) . However lKoenig et al 



(|2008l ) argue that the isolated YSO candidates 
could be genuine YSOs that formed from isolated 
events or that were ejectecj^ from their birthplace 
due to gravitational interactions with other mem- 
bers of their native star cluster. 

We now focus our analysis on the popula- 
tion of clustered YSO candidates since they are 
more reliable tracers of triggered star formation 



^"Assuming a relative velocity of 3 km between a young 
star and its parental cloud, a YSO could travel 6 pc in 
less than 2 Myr. This explains how genuine YSOs can 
contribute to the distributed population in star forming 
regions. 
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events. Various methods can be used to discrim- 
inate between distributed/clustered populations. 
For instance the two-point c orrela tion function 
wa s used by iKarr fc Martini (|2003l ) on W5 and 
by llndebetouw et al.l (12007^ on M16, while the 
minim um spanning t ree t echnique was used on 
W5 bvlKoenig et al.l tOOfj) and on NGC 1333 by 



Gutermuth et al.l (|2008l ). Both methods provide 



valuable qualitative results but dynamical mea- 
surements are absolutely necessary to assign a 
definitive membership to a given YSO. Since such 
data are not available for VulOBl we opt for the 
simplest filtering, namely the nearest neighbor, to 
gain insight into the clustering properties of our 
YSO sample. In practice, we consider all YSO 
candidates and keep only those that have their 
nearest neighbor YSO within 3' (2 pc at 2.3 kpc). 
The result of this selection is a much more clus- 
tered distribution of YSO candidates that traces 
mid-IR morphology even more closely (192 YSO 
candidates are marked as isolated objects). 
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Fig. 6. — Distribution of YSO candidates in VulOBl displayed over the 70 //m negative image from MIPSGAL. Red circles represent 
envelops/ Class I YSO candidates, green horizontal rectangles represents disks/Class II objects, and blue vertical rectangles are optically 
thin disks/Class III objects. Black contours represent the surface density of YSO candidates in VulOBl as described in the text. Their 
distribution is well correlated with the mid-IR morphology. They are found at the tip of most pillars and around Hll regions. The large 
and small dashed circle s around (l.b ) = (59.58,0.12) indicate the diameter of the supernova remnant SNR G59. 5+0.1 as published in 
Tavlor et al. (1992) and Green ( 20061 ) respectively. The orientation and scale of the image are also indicated. 



We also look for trends in the environment of 
YSO candidates compared to field stars. The top 
panel of Figure [7] shows that YSO candidates sit 
preferentially on regions of bright 24 /im back- 
ground emission compared to the more evolved 
stars. Since the 24 ^m emission is mostly due to 
the thermal emission of small dust grains excited 
by UV radiation, the trend for YSO candidates 
to sit on bright mid-IR background implies that 
the identified YSO candidates are largely associ- 
ated with photon dominated regions (PDR) and 
Hll regions. This is consistent with known sce- 
narios of triggered star formation mechanisms and 
with previous stud ies of star forming regions (e.g. 
Koenig et al.l 120081 ). We also compare the visual 
extinction associated with the YSO population to 
the entire population of point sources found in our 
initial catalog (see bottom panel of Figure [7]). We 
find that YSO candidates are mostly located in 
regions of high extinction compared to field stars. 
This is again consistent with the idea that infant 
stars still live in the dense and dusty cloud from 
which they were born. Remarkably the environ- 
mental trends mentioned above are even more pro- 
nounced for the population of clustered YSO can- 
didates (see Figure [7]). 



Note that IChapin et al.l (|2008l ) identified a few 
BLAST sources in the direction of VulOBl that 
are actually located at 8.5 kpc from the Sun (radial 
velocity around —5 km s~^), i.e. in the Perseus 
arm. A dozen of our YSO candidates seem to be 
associated with these distant BLAST sources so 
that these might belong to the Perseus arm. 

5.3. Spectral Energy Distribution 

We have built Spectral Energy Distributions 
(SEDs) for aU of our YSO candidates. Figure M 
presents a small SED sample with the Class des- 
ignation associated with each object. 

We used the SEP fitter tool developed by 
Robitaille et al.l (|2007l ) to extract the physical pa- 
rameters of our YSO candidates. However, the 
limited wavelength coverage (near- to mid-IR) 
resulted in strong degeneracies over the output 
parameters, and we were unable to gain reliable 
information on the physical properties of our YSO 
candidates. 

We looked for possible associations in the 
BLAST field to extend the wavelength coverage 
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Fig. 7. — Histograms of 24 /im background bright- 
ness (upper panel) and visual extinction (lower 
panel) associated with each point sources of our 
initial catalog. YSO candidates appear to sit pref- 
erentially on high 24 /zm background as well as in 
high extinction, i.e. embedded, regions. This is 
even more pronounced for clustered YSO candi- 
dates. (A color version of this figure is available 
in the online journal.) 
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to the sub-millimeter, but in most cases several 
YSO candidates fell within the BLAST beam 
rendering the association ambiguous. For the 
YSO candidates that could be uniquely associ - 
ated with a BLAST source, IChapin et all (|2008l ) 
provide physical parameters derived from SED 
fitting. They find clump masses ranging from 40 
to 400 Mq and temperatures from 19 to 28 K. 
We also use iMolinari et all ( 2008a ) diagnostic di- 
agrams based on MIPS [24-70] and BLAST [250- 
500] color s to identify th e most massive YSOs in 
VulOBl. Molinari et al.l show that the high-mass 
analogues of the low- or intermediate-mass Class 
objects have distinctive MIPS colors. From our 
sample of YSO candidates, we find 6 objects with 
[24 — 70] > 3, which is indicative of an SED peak- 
ing longwards of 70 ^m presumably due to a mas- 
sive infalling envelop. Three of these red MIPS 
sources are located next to the brightest BLAST 
source V30 and are barely discernible on the map, 
two have no BLAST counterparts, and one might 
have a BLAST counterpart but the association is 
ambiguous as another YSO candidate falls within 
the BLAST beam. 

Finally we searched the IRSA for Bolocam 
sources, and we found 24 sources in VulOBl. 
Among them we find isolated sources not asso- 
ciated with any YSO candidates, these could be 
YSOs in their earliest evolutionary phases, i.e. 
starless cores emitting in the millimeter only. The 
other Bolocam sources are usually located close to 
the peaks of the YSO surface density and are asso- 
ciated with groups of YSO candidates. Note that 
almost all Bolocam sources have BLAST counter- 
parts over the Vulpecula BLAST field. 

5.4. Cases of triggered star formation? 

5.4.1. The case of SNR G59. 5+0.1 

Shock waves generated by supernovae can po- 
tentially trigger episodes of star formation at 
the interface between the SNR and the ISM 
200fih . 



(Melioh et al 



The supern ova remnant 



SNR G 59.5-H0.1 was first detected bv lTavlor etld 



()l992l) in the direction of VulOBl at t he posl 



tion (l,h) = (59.58,0.12). iTavlor et al.l describe 
this object as a shell-type SNR with non-thermal 
spectral index {a < —0.4) and a diameter of 15'. 
They also mention the close proximity and pos- 
sible association of the SNR with the Hll region 



Sh2-86. We investigate the possibility of such an 
association and potential signs of triggered star 
formation. 



Guseinov et al.l (j2003l) computed the distance 



to SNR G59. 5-1-0.1 using an empirical formula that 
relates the surface brightness of a SNR (E) with 
its diameter (D). They calibrated the Y, — D re- 
lation against SNRs of known distances, and they 
found that SNR G59.5-f0.1 lies at --11 kpc from 
the Sun. This would imply that SNR G59.5+0.1 
is not associated w ith the Hll regio n Sh2-86. We 
stress however that iGuseinov et al. computed the 
distance to G5 9. 5-1-0.1 based on data from the 
SNR catalog of iGreenI (l2006l) (Y^ighz = 18.1 x 



10^2^ W.m~^.IIz~^.sr~^ and D = 5' from ob ser- 
vations at 1720 MHz), whereas I Taylor et al. re- 
port different values from observations at 327 and 
4850MHz i^iGHz = 0.7x10-21 w.m-2.Hz-i.sr-i 
and D = 15'). We use the VGPS 21-cm data 
to look for traces of the SNR and settle this 
discrepancy. We find a circular structure cen- 
tered at the position of the SNR with a diam- 
eter of 15', which confirms the value found by 



Taylor et al.l . Therefore we compute the di stance 



to SNR G 59.5-F0.1 usi ng iGuseinov etalj 's for- 
mula with iTavlor et al.l 's measurements instead 



of iGreenf s: and we find that the SNR likely lies 
between 2.1 kpc and 5.3 kpc from the Sun pro- 
vided the given uncertainties of the E — _D relation. 
This distance estimate reconciles the hypothesis 
that SNR G59. 5+0.1 is indeed associated with 

S h2-86. 

Reach et al. (I2OO6I) have searched for infrared 



counterparts to the known SNRs in the inner 
galactic plane using IRAC observations from the 
GLIMPSE survey, however they did not detect 
SNR G59. 5+0.1 because of the high confusion 
level in the vicinity of Sh2-86 structured extended 
emission. For the same reason, we could not de- 
tect the SNR on MIPSGAL 24 ^m images either. 
No morphological clues were found in the ancil- 
lary data mentioned in Section 12.21 Neverthe- 
less, the distribution of YSO candidates around 
SNR G59. 5+0.1 reveals two groups of YSOs ap- 
parently lining up along two arcs at the northern 
and southern rims of SNR G59. 5+0.1 (cf Figure|n|) . 
Assuming that the actual diameter of the SNR is 
indeed 15', the two groups do not lie exactly on 
the rim of the shell but slightly outside. We still 
believe that these YSO overdensities could have 
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Fig. 8. — Sample of SEDs from our catalog of YSO candidates. Almost all YSO candidates have detection 
in all four IRAC bands, 75% have detections at 24 iim, 6% at 70 /zm and about half in J, H and K bands. 
Our classification of each YSO candidates is shown on individual plots. The axis scale is the same for each 
subplot. 



arisen from the interaction of the expanding shell 
of SNR G59.5+0. 1 with the neutr al gas of Sh2- 



86. According to IXu et al.l (|2005| ). the age of a 



shell- type SNR such as SNR G59. 5+0.1 with a di- 
ameter of 10 pc (assuming a 15' diameter and a 
distance of 2.3 kpc) ranges from 10'^ to lO'' years 
old. This timescale is comparable to the lifetime of 
the first phases of star formation, making the SNR 
a possible trigger for the surrounding YSO candi- 
dates. These clues tend to confirm the association 
of SNR G59.5-f0.1 with Sh2-86, but since we could 
not find unequivocal evidences for the SNR-HlI re- 
gion association, a fortuitous spatial coincidence 
cannot be excluded. 

5.4-2. Photoevaporation at the tip of the pillars 

An interstellar cloud exposed to the ionizing ra- 
diation of a newly formed star is compressed by 
an ionization-shock front which can focus the neu- 
tral gas into a compact globule. This mechanism 
is called rad iation-driven im plosion (RDI) and is 
described in iBertoldi ( 19891 ). In some cases, this 
leads to the birth of a second generation of stars 
located at the interface of the ionized and neutral 



gas (e.g. Reach et al. 20091 ). 

Even though we could not identify the sculptor 
of the pillars described in Section 13.21 we expect 
the RDI mechanism to be at work in these objects. 
Figure |4] shows the location of our YSO candidates 
in pillars VulPl to VulP6 and VulPll. We find at 
least one YSO candidate at the tip of most of the 
pillars (see Figure |3] and last column of Table . 
YSO candidates can also be found along the pillars 
when the gas is protruding and exposed to the 
ionizing star (e.g. YSO on the West side of VulPS 
on Figure |4]). BLAST also detected sources at 
the tip of every pillar found within the Vulpecula 
BLAST field (e.g. bottom panel of Figure I!]) . 

5.4-3. Embedded star cluster 



Cr404 is a young (9 Myr old, iBica et al.ll2008[ ) 
star cluster located at {l,b) — (59.14,-0.11) em- 
bedded in Sh2-86. We find a high surface den- 
sity of YSO candidates associated with this clus- 
ter. 11 Class n-III and 3 Class I YSO candi- 
dates surround the mid-IR peak emission of Cr404. 
BLAST observations further reveal that the em- 
bedded cluster is coincident with the bright sub- 
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millimeter sou rce V18, also identi fied as IRAS 



19403+2258. IChapin et alJ (|2008| ) classify this 



source as a candidate Hyper Compact (HC) Hll 
region, i.e. the precursor of a high-mass star, and 
they derived a mass of ~150 Mq and a tempera- 
ture of '^SG K for this object. The simultaneous 
presence of a HC Hll region (~ 10^ yr lifetime) 
and Class H YSOs (few 10^ yr lifetime) around a 
young star cl uster such as Cr4 04 is consistent with 
the results of iMolinaril (|2008b[) which suggests that 
the most massive objects in a cluster are the last 
ones to form. 

Other cases of HC Hll regions associated with 
groups of Spitzer-selected YSO candidates are 
found in the BLAST field. For instance, the 
BLAST source V30, the brightest and most m as- 
sive clump according to I Chapin et al. ( 2008 ) is 
coincident with the highest density of YSO can- 
didates in the whole 6.6 square degree map at 
il,b) = (59.80,0.064). 

5.4-4- Large scale propagating Star formation 
Based on radio observations of Hi supershells 



and on energetics arguments, lEhlerova et al 



(|2001[ ) estimate that the Hi shell GS061-h00-h51, 
which encompasses Sh2-87 and Sh2-88, is 6-7 
times larger and 3-4 Myr older than the shell 



GS59.7-0.4-h44 that hosts Sh2-86. lEhlerova et al 
attribute these age and size differences to the de- 
layed expansion of Hll regions, and they suggest 
that star formation might b e propagating from 
Sh2-88 to Sh2-86. Previouslv. fTurneil (|l986[ ) sug- 
gested that the morphology of the three Hll re- 
gions was shaped by a single supernova shock wave 
associated with the fossil Hll region Lynds 792 
{l,b) = (60.81,2.24). Since Lynds 792 is ap- 
proxi mately e quidistant from the three Hll re- 
gions, iTurneil 's scenario implies that star forma- 
tion would be triggered simult aneously, as oppo sed 
to sequentially in the case of lEhlerova et al.l We 
now test the validity of these scenarios using the 
supplemental information gathered in this work. 

Figure [9] shows the evolution of the number of 
YSO candidates as a function of their Galactic 
longitude. The increasing number of YSO candi- 
dates towards Sh2-86 seems to favor the scenario 



of Ehlerova et al.l since the youngest star forming 
region should indeed possess the largest number 
of young stellar objects, assuming identical IMFs 
and instantaneous star formation. Nevertheless 



the shape of this histogram could also be inter- 
preted as Sh2-86 being a more active star forming 
region than the other two Sharpless objects, cre- 
ating more YSOs from a larger molecular reser- 
voir, which is unrelated to any process of propa- 
gating star formation. For a proper interpretation 
of the Figure [9] histogram, we need to compare the 
typical lifetime of a YSO and the timescale neces- 
sary for star formation to propagate in a molecular 
cloud. 



Nomura fc Kamaval (|2001r ) have studied self 



propagating star formation using numerical sim- 
ulations and found that the time delay of se- 
quential star formation sites against the origi- 
nal one. At, correlates with their physical sep- 
aration, Aa:, as At ~ 50M yr\Ax/(0.5kpc)Y'-^ . 
Efremov fc Elmegreen ( 19981 ) have derived a sim- 
ilar expression based on observations of star 
clusters in the Large Magellanic Cloud (At ^ 
26Myr[Ax/{0.5kpc)]°'^). Assuming Sh2-88 and 
Sh2-86 are separated by ^80 pc, we find that star 
formation would take ^1 0-20 Myr to propa gate 
across VulOBl. Besides. lEvans et"aL ( 20091 ) es- 
timate the lifetime of a YSO to be of the order 
of 1-3 Myr based on a large statistical sample of 
YSOs identified as part of the c2d legacy survejf^- 
If star formation was indeed propagating along the 
Galactic plane at the pace derived above, YSOs 
would only be found in a narrow slab of longi- 
tudes since they would age, and fade away in the 
IR, faster than the triggering shock needs to cross 
the molecular complex. We can already rule out 
this scenario based on the presence of YSOs in all 
three Hll regions of VulOBl (see Figure [6]). It is 
possible however that star formation is propagat- 
ing faster than previously estimated due to the 
pronounced inhomogeneities of the propagation 
medium (cf distribution of CO in Figure [1]). If the 
propagation timescale was of the same order of 
the YSOs lifetime, then YSOs would be found all 
across VulOBl with a gradient in the evolution- 
ary phases of the YSO population as a function of 
longitude. 

We use the ratio of Class II-III to Class 0- 
I as an indicator of the aging of the YSO pop- 
ulation. Figure [9] shows that this ratio does 
not present a gradient with respect to longitudes. 
This indicates that if star formation was once 



1 http://irsa.ipac.caltech.edu/data/SPITZER/C2D 
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triggered in VulOBl, then the triggering agent 
would have been independent of Galactic longi- 
tu de, and this def initely rules out the suggestion 
of ihi erova et al.l Nonetheless, the his togram of 
Figure IHI is still consistent with iTurneii 's scenario 
in which star formation occurs as an instantaneous 
burst in the three Hll regions. 

6. Conclusions 

We have presented a thorough description of 
the Vulpecula OB association. We have comple- 
mented the existing observations of isolated ob- 
jects in VulOBl with a global view of the whole 
star forming complex from an infrared perspec- 
tive. We exploited Spitzer legacy surveys MIPS- 
GAL and GLIMPSE data to identify 856 young 
stellar objects with IR-excess emission. We rely 
on the nature and properties of these objects to 
highlight the recent activity of star formation in 
the complex, and we look for evidences of triggered 
star formation. 

We find two populations of YSO candidates: 
one population of distributed objects that likely 
contains IR-bright evolved stars and some gen- 
uine YSOs born in isolation or ejected from their 
parental cloud; and another population of clus- 
tered YSO candidates whose spatial distribution 
correlates very well with the mid-IR morphology 
of the complex. YSO candidates surface den- 
sity peaks locally around the three Sharpless ob- 
jects, the Hyper Compact Hll regions, and other 
embedded star clusters like Cr404. The vigor- 
ous star forming activity around these energetic 
sources is consistent with scenarios of triggered 
star formation mechanisms. Still we cannot as- 
certain that these stars were born as a direct con- 
sequence of their extreme environment, nor that 
they would have formed in the absence of what 
we consider the triggering agents. Nevertheless 
our analysis al l owed to rule out the scenario of 
Ehlerova et al.l (j200l[ ) according to which star for- 
mation was propagating from Sh2-88 to Sh2-86. 
We rather find that the evolutionary stage of the 
YSO population across VulOBl is homog eneous, 
and thus consistent with the scenario of Turner 



(|l98ffl . 

We have further reported the discovery of a 
dozen pillar-like structures in VulOBl, and we 
comment on their morphology from the near-IR 
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Fig. 9. — Top panel: Distribution of class O/I 
and II/III candidates as a function of Galactic 
longitude (bin size ~15'). The number of young 
stellar objects increases towards lower longitudes. 
Bottom panel: Evolution of the class II/III to 
class O/I ratio as a function of longitude. The 
number of class O/I and II/III in each bin is indi- 
cated on the figure. (A color version of this figure 
is available in the online journal.) 
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to the radio regime. We were not able to identify 
the energetic source (s) responsible for the molding 
of the pillars, but we argue that these objects are 
indeed associated with the OB association. Our 
finding of YSO candidates at the tip of most of 
the pillars is consistent with mechanisms of trig- 
gered star formation on small scales. 
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Appendix material 



Table 3 

Pillars identification, coordinates and approximate size. 



Galactic Equatorial 



Pillar ID 1 b RA (J2000) Doc (J2000) Size Position angle 



VulPl 


60. 


.18 


-0.31 


296. 


.36 


23. 


.88 


5-8 


X 


0.5 


250 


VulP2 


60, 


.03 


-0.35 


296, 


.32 


23, 


.74 


1 


X 


0.5 


230 


VulPS 


59, 


.97 


-0.31 


296, 


.25 


23, 


.70 


8 


X 


2 


260 


VulP4 


59, 


.81 


-0.29 


296, 


.15 


23, 


.58 


3 


X 


1 


268 


VulP5 


60, 


.06 


0.10 


295, 


.91 


23, 


.99 


4 


X 


2 


335 


VulP6 


60, 


.11 


0.21 


295, 


.83 


24, 


.10 


5 


X 


2-4 


350 


VulPl 


60, 


.38 


0.21 


295, 


.98 


24, 


.33 


9 


X 


1-4 


17 


VulPS 


60, 


.90 


0.22 


296, 


,26 


24, 


,78 


11 


X 


4 


55 


Vulpg 


61, 


.10 


-0.45 


297, 


.00 


24, 


.61 


4 


X 


1 


158 


VulPlO 


60, 


.77 


-0.63 


297, 


.00 


24, 


.24 


4 


X 


2 


190 


VulPll 


60, 


,21 


-0.45 


296, 


,52 


23, 


.84 


4 


X 


2 


280 


VulPl 2 


59, 


,60 


-0.14 


295, 


,90 


23, 


.47 


3 


X 


1 


133 


VulPl 3 


59, 


,51 


-0.22 


295, 


,92 


23, 


.36 


3 


X 


0.5 


160 


VulPli 


59, 


.49 


-0.19 


295, 


.88 


23, 


.35 


2 


X 


0.5 


150 



Note. — Coordinates and position angles are in units of degrees. Pillar sizes are in units of 
arcminutes. Position angles are given between the long axis of the pillar and the West axis, 
positive values are towards the North. 
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Table 4 

Morphological description of the pillar structures identified in VulOBI. 





MIPS 






IRAC 






Pillar ID 


[70] 


[24] 


[8.0] 


[5.8] 


[4.5] 


[3.6] 


YSO'' 


VulPl 


Bright Pillar 
Compact Source 


Structured 


Structured 


Faint 


None 


Faint 


Y 


VulPZ 


Faint 


Faint 


Faint 


None 


None 


None 


N 


VulPS 


Bright Pillar 
Compact Source 


Bright 
Structured 


Bright 
Structured 


Faint 


Faint 


Faint 


Y 


VulP4 


Faint Tip 


Faint 


Faint 


None 


None 


None 


Y 


VulPS 


Faint Pillar 


Bright Pillar 


Bright Pillar 


No Pillar 


No Pillar 


No Pillar 


Y 


Pedestal 


Pedestal 


Pedestal 


Faint Pedestal 


Faint Pedestal 


Faint Pedestal 






ramT; 


ramx 


raint 


None 


None 


None 


V 


VulP7 


Faint 
Structured 


Bright 
Structured 


Faint 
Structured 


None 


None 


None 


Y 


VulPS 


Faint Tip 


Faint 
Structured 


Faint 
Structured 


None 


None 


None 


Y 


VulP9 


Faint 


Bright tip 


Bow Shock 


Faint Tip 


None 


Faint 


N 


VulPlO 


Faint 


Bright Tip 
Possible Jets 


Faint 


None 


None 


None 


Y 


VulPll 


Bright Tip 


Bright Tip 


Bright Tip 


Faint Tip 


Faint Tip 


Faint Tip 


N 


VulPl 2 


Bright Tip 
Pedestal 


Bright Tip 
Pedestal 


Bright Tip 
Pedestal 


Faint 


Faint 


Faint 


Y 


VulPl 3 


Faint 


Bright Edges 


Bright Edges 


None 


None 


None 


Y 


VulPU 


Bright Tip 


Bright Tip 


Bow Shock 
Bright Tip 


Bow Shock 


Bow Shock 


Bow Shock 


Y 



^Indicates the presence of a YSO candidate at the tip of the pillaj:. 

Table 5 

List of point sources excluded from the YSO catalog and identified as probable 
contaminants based on their IRAC fluxes. 



Galactic 2MASS IRAC MIPS 



GLIMPSE source name 


Glon 


Glat 


J 


H 


Ks 


[3.6] 


[4.5] 


[5.8] 


[8.0] 




[24] 


[70] 


Type 


SSTGLMC G059.6291-00.8301 


59.6291 


-0.8301 


12.97 


12.59 


12.35 


11.96 


11.63 


11.41 


9.84 








Galaxy 


SSTGLMC G059.7319+00.2144 


59.7320 


0.2144 


12.95 


12.47 


12.23 


11.98 


11.71 


11.68 


10.23 


7 


77 




Galaxy 


SSTGLMC G059. 9431+00. 1783 


59.9432 


0.1783 




15.06 


14.19 


13.02 


12.40 


11.68 


9.76 


6 


19 




Galaxy 


SSTGLMC G060. 9219-00. 1041 


60.9220 


-0.1041 


14.37 


13.44 


12.98 


12.15 


11.87 


11.83 


10.53 








Galaxy 


SSTGLMC G061. 1282-00.7549 


61.1282 


-0.7549 


13.60 


12.59 


12.32 


12.11 


12.14 


11.77 


10.63 








Galaxy 


SSTGLMC G061. 645.3-00. 7439 


61.6454 


-0.7439 








12.19 


11.91 


11.79 


9.53 








Galaxy 


SSTGLMC G060. 0154-1-00. 1111 


60.0154 


0.1111 








12.88 


10.83 


10.01 


9.56 






-4.96 


Shock 


SSTGLMC G059.6940-I-00.1840 


59.6940 


0.1840 








12.80 


11.62 


11.14 


10.69 






-4.82 


Shock 


SSTGLMC G059.7893-I-00.6298 


59.7893 


0.6298 








12.61 


10.75 


10.19 


9.27 






-4.81 


Shock 


SSTGLMC G059.3780-00.2438 
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Fig. 10. — Composite image of the Hll region Sh88B as seen by Spitzer (blue, green and red are [3.6], [5.8] 
and [8.0] respectively). 
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Fig. 11. — Composite image of the Hll region Sh86 as seen by Spitzer (blue, green and red are [5.8], [8.0] and [24] respectively). 
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